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CHAPTER  I 


Introduction 

A  maximum  flow  network  it  defined  by  a  eet  of  arc*  and  a  eet 
of  points  called  nodes.  Each  arc  joins  two  nodes  and  has  associated 
with  it  a  positive  capacity  which  represents  the  maximum  amount  of 
flow  that  may  pass  over  it.  One  of  the  nodes  is  designated  as  the 
source  and  another  as  the  sink.  From  these  nodes,  arc,  and  capaci¬ 
ties  the  maximum  amount  of  flow  that  may  pass  from  source  to  sink 
may  be  calculated. 

This  investigation  is  concerned  with  a  sensitivity  analysis  on 
such  networks.  Specifically,  each  arc  of  the  network  is  assumed  to 
be  subject  to  breakdowns  which  result  in  a  reduction  in  its  capacity. 
The  problem  is  to  find  the  greatest  reduction  in  maximum  How  pos¬ 
sible  if  n  breakdowns  occur  and  to  find  where  these  breakdowns 
must  occur  for  this  reduction  to  result.  The  methods  developed  in 
this  report  solve  this  problem  for  a  certain  class  of  networks  known 

as  planar  networks. 

The  first  method  solves  the  problem  exactly  where  it  is  as¬ 
sumed  that  each  arc  is  subject  to  only  one  breakdown  and  the 
amount  by  which  the  capacity  of  an  arc  is  reduced  due  to  a  break¬ 
down  is  a  deterministic  quantity.  The  algorithm  developed  could, 
however,  easily  be  modified  to  aUow  for  multiple  breakdowns  on  one 

arc. 

The  second  method  solves  this  problem  approximately  when 
the  amount  by  which  the  capacity  of  an  arc  is  reduced  is  a  random 
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variable  with  unknown  distribution  but  with  known  mean  and 
variance.  For  this  case,  multiple  breakdowns  are  allowed  for  the 
individual  arcs.  The  algorithm  developed  can  also  be  used  to  solve 
the  problem  more  exactly  when  these  random  variables  have  normal 
distributions. 
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CHAPTER  II 


Basic  Concepts 

1.  Max-Flow  Min-Cut  Theorem 

The  central  idea  in  the  solving  of  maximum  flow  network 
problems  is  summarized  in  the  max-flow  min-cut  theorem 

Definition:  Consider  a  network  consisting  of  nodes  which  in¬ 
clude  a  source  and  a  sink,  and  capacitated  arcs  which  join  two  nodes. 
Let  A  and  B  be  a  partition  of  the  nodes  such  that  the  source  is 
in  A  and  the  sink  is  in  B  .  Then  the  set  of  arcs  which  join  a 
node  in  A  to  a  node  in  B  is  called  a  cut  set  and  is  denoted 
[A,  B]  Furthermore,  the  value  of  this  cut  set,  V[A,  B]  ,  is  e- 
qual  to  the  sum  of  the  capacities  of  its  arcs. 

A  property  of  any  cut  set,  [A,  B]  ,  is  that  any  path  from 
source  to  sink  must  use  at  least  one  of  its  arcs.  Thus,  it  is  appar¬ 
ent  that  the  maximum  flow  cannot  exceed  the  minimum  value  of  all 
cut  sets.  The  max-flow  min-cut  theorem  states  that  the  maximum 
flow  actually  equals  the  minimum  cut.  ^ 

2.  The  Topological  Dual 

The  topological  dual  of  a  network,  when  defined,  is  another 
network  in  which  the  arcs,  instead  of  having  capacities,  have 
lengths.  Furthermore,  there  is  a  one-to-one  correspondence  be¬ 
tween  the  proper  cuts  of  the  original  network  and  the  routes  through 
the  dual,  and  the  problem  of  finding  the  minimum  cut  may  be  re- 

^Reference  8. 
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duced  to  one  of  finding  a  shortest  route. 

Let  the  original  maximum  flow  network  be  called  the  primal. 
To  the  primal  add  an  artificial  arc  extending  from  the  source  to  the 
sink  and  having  a  capacity  of  zero.  The  resulting  network  will  be 
referred  to  as  the  modified  primal.  The  dual  network  is  defined  if 
and  only  if  the  primal  is  source- sink  planar ,  a  source- sink  planar 
network  being  one  where  the  modified  primal  can  be  drawn  on  a 
^here  in  such  a  way  that  no  two  arcs  intersect  except  at  a  node. 

When  defined  the  dual  is  constructed  in  the  following  manner: 

1.  Draw  the  modified  primal  on  a  sphere  in  such  a  way  that  no 
two  arcs  intersect  except  at  a  node. 

2.  Place  a  node  in  each  mesh  of  the  modified  primal.  Let  the 
node  in  one  of  the  two  meshes  bounded  by  the  artificial  arc  be  the 
source  and  the  node  in  the  other  of  these  two  meshes  be  the  sink. 

3.  For  each  arc  except  the  artificial  one  construct  an  arc  of 
the  dual  that  intersects  it  and  joins  the  nodes  in  the  meshes  on  either 
side  of  it. 

4.  Assign  each  arc  of  the  dual  a  length  equal  to  the  capacity 
of  the  primal  arc  it  intersects. 

An  example  of  a  network  and  its  dual  is  shown  in  Figure  1. 
Letting  a  route  through  the  dual  be  any  path  from  its  source  to  its 
sink,  it  follows  that  there  is  a  one-to-one  correspondence  between 
the  proper  cuts  of  the  primal  and  the  routes  of  the  dual.  ^  Specifi¬ 
cally  if  A  is  any  route  through  the  dual  then  the  arcs  of  the  primal 

TIT - 

1  'Reference  8. 
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which  intersect  A  form  a  proper  cut  and  conversely  if  B  is  a 
proper  cut  of  the  primal,  then  the  arcs  of  the  dual  which  intersect 
B  form  a  route. 
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CHAPTER  in 


Deterministic  Case 

1.  Problem  Formulation 

The  problem  solved  in  this  chapter  is  the  following: 

A  source- sink  planar  maximum  flow  network  is  given  in  which 
each  arc  is  subject  to,  at  most,  one  breakdown:  a  breakdown  re¬ 
sulting  in  a  reduction  in  the  capacity  of  that  arc  by  a  known  quantity. 
It  is  desired  to  find  the  greatest  reduction  in  maximum  flow  possible 
from,  at  most,  n  breakdowns  and  where  these  breakdowns  must 
occur,  to  give  this  reduction. 

Since  the  maximum  flow  of  a  network  is  equal  to  its  minimum 
cut  which  in  turn  is  equal  to  the  length  of  the  shortest  route  through 
its  dual,  it  follows  that  solving  this  problem  is  equivalent  to  solving 
the  following  problem  for  the  dual: 

A  network  is  given  in  which  the  arcs  are  assigned  lengths. 
Each  arc  is  subject  to  an  improvement  which  results  in  a  decrease 
in  its  length  by  a  known  quantity.  It  is  desired  to  find  the  smallest 
value  possible  for  the  length  of  the  shortest  route,  if  no  more  than 
n  improvements  occur  and  where  these  improvements  must  occur 
in  order  to  achieve  this. 

It  is  the  latter  problem  that  the  algorithm  of  this  chapter 
solves.  Throughout  the  rest  of  this  chapter  all  nodes  and  arcs  re¬ 
fer  to  the  dual  network  unless  specified  otherwise. 
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2.  The  Algorithm 

Let  an  i-arc  path  to  node  a  be  any  path  from  the  source  to 

node  a  along  which  at  most  i  improvements  occur.  Also,  let 

L  (to  be  determined  later)  be  the  length  of  the  shortest  i-arc 
a,  i 

path  to  node  a  .  The  algorithm  assigns  to  each  node  a  ,  n  +  1 
labels,  (D,t,k)a  Q  ,  (D,t,k)a  x  .  ,  (D,t,k)a  n  .  The  individ¬ 

ual  components  of  (D,  t,k)a  ^  will  be  denoted  by  Da  .  ,  ta^ 

and  k  .  .  Initially,  it  is  known  that  all  D  .  >  L  .  The 

a,  l  a,  i  —  a,  i  -* 

O  are  then  decreased  in  such  a  way  as  to  preserve  this  initial 
a,  i 

property  until  all  Da  .  =  La  .  '  .  At  this  point  Dg  n  is  the 

length  of  the  desired  path  where  5  is  the  sink.  The  components 

t  and  k  are  tracers  which  are  used  to  find  the  desired  path 
a,  i  a, l 

itself. 

Let  l(a,b)  be  the  length  of  arc  (a,b)  ,  d(a,b)  the  decrease 
in  this  length  due  to  an  improvement,  S  the  source,  and  5  the 
sink.  Note  that  l(a,h)  =  l(b,a)  and  d(a,b)  =  d(b,a)  The 
algorithm  for  finding  the  length  of  the  desired  route  is  as  follows: 

1.  For  i  =  0,1,  .  . .  ,n  set  Dc  .  =  0  and  D  .  =  co  for 
a  £  S  .  Set  m  =  0  . 


2.  Check  each  arc  (a,b)  and: 
a.  If  D  >  D.  +  l(a,b) 

& i  XXX  Df  XXX 

D»  «  =  Dh  m  +  M*.b) 
a,  m  o,  m 

t  =  b 

a,  m 

k  =0 

a,  m 


set 
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b.  If  m  >  3  and  D  >0,  ,  +  I(a,  b)  -  d(a,b) 

—  a,  m  d,  m- 1 

<•  ^  ■ 

w  \-r  4- 

set 

D  _  =  D,  ,  +  l(a,b)  -  d(a,  b) 
a ,  m  o ,  m-i 

t  =  b 

a,  m 

k  =  1 

a,  m 

3.  Repeat  ~2  until  no  more  changes  can  be  made.  Then  if 
m  <  n  ,  increase  m  by  1  and  go  back  to  2  If  m  =  n  , 
terminate,  as  Dg  n  is  the  length  of  the  desired  route. 

The  desired  route  itself  may  be  found  by  the  following  proce¬ 


dure: 


1.  Set  m  =  1  ,  =  a 


1 


li =  n 


2.  Let  am+1  =  t  and  W  =  V  -  i  ■ 

mm  mm 

3.  If  am+j  £  S  ,  increase  m  by  1  and  go  back  to  2  . 
Otherwise  terminate. 


S  =  a 


. a^  =  3  is  the  desired  path.  The  improvements 


must  occur  on  those  arcs  (a.  ,,  a.)  where  k  .  =  1  The  arcs 

JTl  J  aj,  lj  * 

in  the  primal  on  which  the  breakdowns  must  occur  are  the  arcs 
which  inter sectlthese  arcs  of  the  dual. 


3.  Justification  of  the  Algorithm 

The  procedure  for  justi^ring  the  algorithm  of  the  last  section 
will  be  to  show  first  that  all  D  .  >  L  .  always  and  D  .  =  L  . 

1  “  A)  1  l  •  i  • 

at  termination.  Then  a  relationship  between  the  L  .  will  be 

1 

established  to  help  verify  the  process  of  tracing  out  the  desired  path. 
Lemma  1:  D  .  >  L_  .  all  a  and  all  i  . 


y 
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Proof:  Assume  that  at  one  point  D  .  >  L  .  all  a  and 
-  ‘  a,  i  —  a.  i 

all  i  .  Suppose  that  some  .  is  changed  to  Da  ^  •  Then 

there  is  a  node  b  such  that  either: 

1.  D  .  s  Dv  .  +  l(a,b)  or 

21)  1  D,  1 

2.  0  .  =*  D.  .  ,  +  l(a, b)  -  d(a, b)  . 

a,i  b, l-l 

In  the  first  case  any  i-arc  path  to  node  b  of  length  L^  ^  plus 

arc  (a,  b)  at  its  true  length  of  l(a,b)  is  an  i-arc  path  of  length 

less  than  or  equal  to  D  In  the  second  case  any  (i  -  l)~arc 

a,  l 

path  to  node  b  of  length  L^  .  ^  plus  arc  (a,  b)  at  its  improve 

ment  length  of  l(a,b)  -  d(a,b)  is  an  i-arc  path  to  node  a  of 

length  less  than  or  equal  to  .  .  Thus  the  relationship 

D  .  >  L  all  a  and  all  i  still  holds.  Initially, 
a,  l  —  a,  l 

D0  .  =  0  -  Le  .  and  D  .  =  oo  >  L  .  for  a  *  S  .  The  lemma 
S.i  S,  i  a,i  a,  l 

follows  from  induction. 

Lemma  2:  After  a  finite  number  of  iterations,  Da  q  =  q 

for  all  a  and  remains  at  that  value  for  all  subsequent  iterations. 

Proof:  Let  S,a. . a  ,  a  be  a  shortest  route  from  S 

. ■■■  i  m 

to  a  .  After  one  examination  of  the  arcs, 


Da.,0  -  1<s-al) 
After  2  iterations, 


1(S,  a^)  +  1  (a^,  a^) 


After 


m  +  1 


iterations, 

Da,  q  <  1(S,  a^) 


m-1 

+  £  1(ai-  ai+l>  +  1(am'  a)  =  La,  0 

i=l 
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Of  course  D  n  will  remain  at  L  ^  tor  all  subsequent  itera- 

a  t  u  a  i  w 

tions  since  the  ^  are  non- me  leasing  and  q  >  Lft  q 

lemma  >.  Thus  the  lemma  holds  for  any  particular  node.  Let 

m(a)  be  the  number  of  iterations  required  in  order  that  q  =  J-<a  q 

Then  after  max  m(a)  iterations  D  ~  =  L  .  for  all  a  . 

a  a ,  u  a  <> 


Lemma  3:  Suppose  5,3,,  ..,a  is  a  shortest  i-arc  path 

— - — —  i  m 

from  S  to  a  If  arc  (a  ,,a  )  has  an  improvement  in 

m  m-i.ni 


this  path  then 
1.  L 


.  =  L  .  ,  +  if  a  -.a  )  -  d'a  a  ) 

a  ,i  a  ,,i-i  m- m  m-1  m 

m  m-i 


Otherwise 


2.  L  =  L 
a  ,  i 
m 


+  l(a  ..a  I 
a  , ,  i  m-1  m 
m-1 


Proof:  Suppose  arc  (a  ,fa  )  has  an  improvement  in  this 
-  m-  i  m 

path.  Then  S.a. . a  ,  is  an  (i-U  -arc  path  to  a_  .  of 

r  1  ml  r  m-I 

length  L  -  Ha  ,  ,a  )  +  d{a  ,,a  )  and 

a  ,  i  m-i  m  m-i  m 

m 

L  <.  L  .  -  Ua  ,,a  >  +  d'a_  ,,a  )  or 

a  ,,i  -  a  ,i  m-i  m  m-1  m 

m-1  m 

L  .  >  L  .  ■+  1(4  ii  4  )  -  d(a 
m  m-1 

Furthermcre,  any  (i-l)-arc  path  of  length  L 


to  a 


•  t  vv  c a  « 

a  ,,i-l  m-1 

m-1 


and  arc  fa^  at  lts  improvement  length  is  an  i-arc  path 

to  a^  of  length  .  ,  +  l(a^  .,a  )-d(a  ,,a  )  and  the 

m  a  , ,  l- .  m-i  m  m-i  m 

m-1 

first  equalitv  holds  If  arc  (a  ,,a  )  does  not  have  an  improve- 

m  -  i  m 

ment,  then  S,a,,.  .,a  .  is  an  i*arc  path  to  a  ,  of  length 

1  m-1  m-i 


-  il  - 


L  .  -  l(a  ,,a  )  and  L  .  >  L  .  +  l(a  ,,a  )  . 

a  ,1  '  m-1  m  ,  1  —  a_  ,,1  m-1  m 

m  m  m-l 


and  arc 


Also,  any  i-arc  path  to  a  ,  of  length  L  . 

m-1 

(a  ,,a  )  is  an  i-arc  path  of  length  L  .  +  l(a  ,,a  ) 

m-l  m  a _  «» i  m-i  m 

m-i 

the  second  equality  holds. 


and 


Lemma  4:  Suppose  S,a. . a  is  the  shortest  i-arc  path 

■  ■■  — l  m 

from  S  to  a  Let  (a  ,a  .)  be  the  last  arc  in  this  path  to 

m  r  r+i 

have  an  improvement.  Then,  if  i  >  1  : 


m-1 


j  .  =  L  .  ,  -  d(a  , a  ..)  + 
a  ,i  a  ,i-l  '  r  r+r 
m  r 


1  l(aj-Vi) 


Proof: 


-  1 


V*  •  V  *  . . . 

La  ..i  '  La  ,  i-1  =  1(Var+l'  *  d(ar,ar+l’ 
r+1  r 

from  Lemma  3.  Summing  these  equations,  gives 

m-1 

L  .  =  L  .  ,  -  d(a  ,a  L.)  +  /  l(a.,a.  ,) 

a  ,i  a  ,  i-I  '  r  r+i  /  j  i+I 
m  r  ^  j  j 

J=r 

Theorem  5:  The  algorithm  terminates  after  a  finite  number  of 

iterations  with  D  .  =  L  all  a  and  all  i 

a,  i  a, l 

Proof:  Suppose  that  after  a  finite  number  of  iterations 

O  .  =  L  all  a  and  all  i  <  M  .  Let  S,a,,  .  .  .  ,a  be  an 
a,  i  a,  i  —  1  m 

(M+l)-arc  path  of  length  Lft  '■  and  let  (ar*ar+i)  t*ie 

last  arc  in  this  path  in  which  an  improvement  occurs.  Then,  after 
one  additional  iteration, 

Da  rM+l  -  La  ,  M  d*ar,ar+l*  +  1(ar'ar+l* 
r+1  r 
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and  after  m  -  r  additional  iterations, 


m-l 

Da  ,M+1  -  La  .M  '  d(ar>ar+l)  +  £  1(aj,aj+l)  =  La  . 
m  r  j  1X4 


M+l 


and  the  theorem  holds  for  any  particular  .  Let  m(a) 

be  the  number  of  iterations  required  in  order  that  Da  = 

La.  M+l  •  The"  aft*r  mix  m(a)  iterationa  Da.  M+l  *  La.  M+l 

all  a  .  Since  D  „  =  L  „  after  a  finite  nSmber  of  iterations, 

a ,  0  a,  u 

the  theorem  follows  from  induction. 

Of  course,  any  route  from  source  to  sink  with  n  or  fewer 

improvements  is  an  n-arc  path  to  5  .  Thus  at  termination 

W  is  the  length  of  the  desired  route.  It  now  remains 
S i  n  Sj  n 

to  justify  the  procedure  for  tracing  out  the  desired  path  itself. 


Lemma  6:  If  k  .  =  1  ,  then 

j’j 


*  •  -  L  .  +  l(a.,a.  .)  -  d(a.,a.  .) 

j,lj  aj+i,lj+i  1  J  3  J 


and  if  k  =  0  then 
j  j 

aj’lj  Vi'Vi  J  J+l 

Proof:  Suppose  k  .  =  1  . 

J  J 

Then 

L  —  D  .  +  l(a.,a.  .)  -  d(a.,a.  .) 

aj,lj  aj+r  j+i  J  J+1  J  J+1 

held  when  the  label  L  was  assigned.  Furthermore, 

d  •  i  1 1 

J  J 

D  =  L  at  this  time  for  if  not  any  i.  ,  =  (i. 

V-Vi  aj+r  lj+i  Jfl  J 


-  l)-arc 
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path  toAiode  a.,,  of  length  L  .  plus  arc  (a., a.  ,)  at 

j+i  aj+rlj+i  J  J+1 

its  improvement  length  is  an  ij  -  arc  path  to  node  aj  of  length  less 
than  L  .  .  Also,  if  k  .  =  0  »  then 

Vj  Vj 

L  i  =  D  .  +  l(a.,a.  j) 

j*  j  j+r  j+1  }  J 

held  when  the  label  L  was  assigned.  Furthermore 

•ii  ij 


j-  ‘j 


s  L_ 


.  -  *4  .  held  for  if  not  any  i.  ,  =  i.  -  arc  path  to 

‘j+r  *j+i  *j+r  j+i  J'''1  J 

node  a^+j  plus  arc  (a^.a^)  at  its  true  length  of  l(aj,a^|)  is 

an  i.  -  arc  path  to  node  a.  of  length  less  than  L  . 

j  j  *j.*j 

Lemma  7:  The  procedure  for  tracing  the  desired  path  is 

finite,  (i.  e. ,  there  exists  r  such  that  i(r)  *  S  '#) 

Proof:  L_  *  <  L_  a  and  i4^,  <  i4 
—  j+i*  j+i —  j*  j  J  1 

with  strict  inequality  for  at  least  one  of  these.  Thus  no  pair  (ar,  iy) 
can  be  repeated  and  the  procedure  must  be  finite. 


Theorem  8:  Let  a^, . . . , 
tracing  procedure.  Then  S*i 
ment  occurring  on  arc  (a^a^ 


aa^|  be  the  nodes  found  by  the 
. . . ,  a^  «  S  with  an  improve- 
)  if  and  only  if  k  .  =  1  is  the 


shortest  n-arc  path  to  5. 


Proei:  Suppose  S  -  »m+i; ; ; .  f. 

arc  (a., a..,)  if  and  only  if  ka  .  »  l 
J  J+l 

Then: 


of  length  La  . 

V  *» 


with  an  improvement  on 
is  an  if  -  arc  path  to  ay 
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1.  L  =  L  .  +  l(a  ,,  a  )  if  k  . 

a  iii  *  a  ,1  r-l  r  a  lti  « 

r-1  r-1  r’  r  r-1  r-1 


=  0 


or 


2-  K  i  =L»  5  +  l(a  ,.a >  -  d(a  ,,a  )  if 

ar-l’  lr-l  Vlr  r  1  r  r  1  r 


k  =  1 

ar-l’lr-l 


In  either  case  a 


m+r  •  *  •  *ar  1  an  improvement  on  arc 


(a^.a^)  if  and  only  if  ka  .  =  1  is  an  i^^-arc  path  to  ar.j 

J  J 


.  Since  a  ,,  =  S  is  a  path  of  zero  length 
ar-l’  lr-l  m+1 


of  length  L 

and  hence  an  i  ..-arc  path  to  a  xl  of  length  L  . 

m+1  m+1  am+l’  lm+l 

it  follows  from  induction  that  am+i*  •  •  •  *  aj  is  an  n-^rc  path  to 

a.  of  length  L  .  =  L*  . 

1  ®  »!•»!  s-n 

This  completes  the  justification  of  the  procedure  for  finding 

the  desired  route  itself.  The  arcs  where  improvements  must  occur 

are  those  arc  8  (a.,  a.,,)  where  k  .  =1  .  The  arcs  of  the 

j  j+i  «j-ij 

primal  upon  which  breakdowns  must  occur  in  order  to  minimize 
the  maximum  flow  are  those  which  intersect  thesefarcs  of  the  dual. 
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CHAPTER  IV 


Stochastic  Case 

1.  Problem  Formulation 

The  arcs  of  the  maximum  flow  network  in  thi6  chapter  are 
subject  to  anywhere  from  0  to  n  breakdowns.  The  decrease 
in  capacity  of  arc  (a,b)  due  to  i  breakdowns  is  a  random  variable 
with  unknown  distribution  but  with  known  mean  and  variance.  The 
mean  and  variance  of  this  quantity  will  be  denoted  as  p(a,b,i) 
and  <r^(a,b,i)  respectively.  Furthermore,  these  distributions 
are  independent  for  the  different  arcs.  It  is  desired  to  find  the 
smallest  possible  value  of  F  resulting  from,  at  most,  n  break 
downs  where  F  satisfies:  % 

P  {max  flow  >  F}  <  /3 

where  P  stands  for  probability  and  is  a  small  number  be¬ 

tween  zero  and  one. 

This  problem  can  also  be  formulated  in  terms  of  the  dual  net¬ 
work.  Specifically,  each  arc  of  the  dual  is  subject  to  anywhere  from 
0  to  n  improvements.  The  decrease  in  length  of  arc  (a,b)  due 
to  i  improvements  is  a  random  variable  with  mean  p(a,b,i)  and 
variance  <r*(a,b,  i)  It  is  desired  to  find  the  smallest  possible 
value  of  F  resulting  from,  at  most,  n  total  improvements  on 
n  or  less  arcs  where  F  satisfies: 

P  {min  route  >  F)  <  /S 

It  is  the  latter  formulation  that  is  considered  in  this  chapter. 
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Since  examples  can  be  constructed  to  show  that  different  distribu- 
tions  with  the  same  |±(a,b,i)  and  <r2(a,b,i)  can  have  different 
F  and  different  locations  for  the  improvements,  the  information 

l 

given  is  not  sufficient  to  determine  a  solution.  Accordingly,  the 
distribution  considered  will  be  the  one  whose  smallest  value  of  F 
is  maximum. 

Suppose  that  if  improvements  occur  on  certain  arcs,  the  length 

of  a  particular  route  is  a  random  variable  L  with  mean  p  and 
2 

variance  <r  It  follows  from  Tchebysheff  s  extended  lemma^ 


P  {L  -  p  >  e}  <  * 


-  JZ"  '2 

€  +  (T 


or 


P  {L  >  p  +  e}  < 


2 

a 

;z~  2 

e  +  <r 


If  it  is  required  to  minimize  p  +  e  under  the  condition  that: 

2 


;z-  2  - 

e  +  <r 


<  P 


it  follows  that 


and  one  obtains 


P  {L  >  p  +  -  1  }  <  0 

Thus  the  quantity  p  +  tryi  -  1  play,  the  role  of  the  length 


in 


Reference  9,  pp.  111-126. 
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of  the  above  mentioned  route  and  will  be  referred  to  as  its  effective 
length.  The  algorithm  to  be  presented  in  this  chapter  finds  those 
arcs  upon  which  the  n  improvements  must  occur  in  order  that  the 
minimum  effective  length  of  all  routes  through  the  dual  be  minimized. 

In  order  that  this  be  a  valid  criteria  the  following  assumptions 
are  made: 

1.  p(a,b,i)  <  l(a,b)  and  is  strictly  increasing  in  i  . 

2.  p(a,b,i)  -  <r(a,  b,  i)"^  “  1  >  0  for  i  >  1  . 

3.  p(a,b,i)  >  p(a,l>,T)  implies 

p(a,b,  i)  -  <r(a,b,i)  Yp  -  1  >  p(a,T>,T)  -  <r(a,T>.T)  x^p  -  1 

These  assumptions  assure  that  the  solution  to  this  problem  involves 
exactly  n  breakdowns,  assure  that  the  effective  length  of  an  arc 
is  always  greater  than  zero  and  generate  arguments  in  favor  of  the 
efficiency  of  the  algorithm. 

2 .  The  Algorithm 

Let  an  i-arc  path  to  node  a  be  any  path  from  the  source  to 

node  a  with  i  or  fewer  improvements  occurring  on  its  arcs. 

Eacl$node  is  assigned  n  +  1  sets  of  labels,  the  set  numbers  being 

designated  as  0, 1,  ....  n  .  Each  label  consists  of  four  components 

and  is  denoted  (p,<r  ,  t,k)^  .  where  a  designates  the  node,  i 

a,  i 

the  set  number  and  j  the  rank  within  the  set.  The  individual  com- 
ponents  are  denoted  as  M-i  •  »  i  »  an<*  k~  .  .  The 

alyl  CL  I  1  cl  9  1  CL  9  1 

*  O  ' 

quantities  p^  .  and  (<r  are  related  to  the  mean  and  vari- 

a,  i  a , i 

ance  of  i-arc  paths  to  node  a  and  at  termination  pg  n  and 
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« 


7  *r 

(o“yJ-  ^  are  the  mean  and  variance  of  the  length  of  the  desired  path 
where  7  =  max  j  in  this  set.  The  components  .  and  . 
are  tracers  which  are  used  to  find  the  path  itself.  Again,  S  and 
3  refer  to  the  source  and  sink  respectively. 

The  algorithm  for  finding  the  mean  and  variance  of  the  length 
of  the  desired  route  is  as  follows: 


for 


1.  Set  i  =  0 

2.  Set  ps  ^  )S  i  =  0  .  Ka>i 
a  £  S  and  s  =  i  . 


/  2,1 
(<r  ) 
x  'a,  l 


=  00 


3.  Consider  each  a  and  all  arcs  of  the  form  (a,b)  . 


For  each  j  consider  the  quantity: 


L  s  .-8  +  l(a,b)  -  p(a,b,  s) 

,b* 


and  delete  from  the  set  i  at  node  a  all  labels  j'  which 


satisfy: 


+  l(a,b)  -  (i(a,b,  *) 

/ .  +  P’ .  >  L 

^a, i  a, i  VP  - 


with  strict  inequality  for  at  least  one  of  these.  If  there  is  no  j1 
satisfying: 

*  9 

i  -  ^b.i-s  +  l<a’b)  *  ^a*b»  ®) 

pj'  .  +  <rj'  .  l/i  -  1  <  L 
*a,i  a,i  VP 

introduce  the  label  (p,<r^,t,k)  .  where 

a,  l 
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|t"  .  =  +  Ma.b)  -  |x(a,b,  a) 

I'2)!. *Z(a.b.B) 


If  any  changes  result  in  the  labels  after  examining  an  arc,  rej-rank 
the  labels  in  order  of  increasing  first  component.  This  results  in  a 
change  in  the  superscript  j  for  some  labels.  In  this  set,  j  now 
takes  on  the  values  1,  2,  .  .  . ,  r  where  r  is  the  number  of  labels 
in  set  i  at  node  a  . 

4.  If  s=0  go  to  5  .  Otherwise,  decrease  s  by  1 

and  go  back  to  3  . 

5.  Repeat  3  until  no  changes  in  the  labels  of  set  i  at 
node  a  result. 

6.  If  i  <  n  ,  increase  i  by  1  and  go  back  to  2  If 

i  s  n  ,  terminate,  as  n  and  (or  ^  are  the  mean  and  vari¬ 
ance  of  the  length  of  the  desired  path  where  J  equals  the  number 
of  members  of  set  n  at  3  . 


The  path  itself  may  be  found  by  the  following  procedure: 

1 .  Set  m  =  1  ,  a  j  =  *5  ,  and  =  n  . 

2.  Let  (|i,  <r  ,  t,  k)^ry  ^  be  the  label  of  highest  rank  among 


set  n  at  5  . 


3.  Set  k  =  k  m  . 

m  a  ,  l 
m  m 


4.  Increase  m  by  1  .  Set  i _  =  i _  ,  -  k  .  and 

m  m-i  m-i 
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i 


a  =  t-1"1-1  . 
m  a  ,,  1  , 

m-1  m-1 


If  a  =  S  ,  go  to  6  .  Otherwise,  find  a 
m 


label  in  set  i  at  node  a  which  satisfies^ 

m  m 


i  at  uuuc  a  w  uui  baubucb 

m  m 

«  • 

pJm  .  =  fi^’1  -  l(a  ,,a  )+  |x(a  ,,a  ,k  .) 

ra  ,  i  a  ,,i  ,  '  m-1  m  m-1  m  m-1 

m  m  m-1  m-1 

(<rVm  .  =  -  <r^( a  ,,a  ,k  .) 

'  a  ,i  '  a  ,  '  m-1  m  m-1 

m  m  m-1  m-1 

5.  Go  back  to  3 

6.  Terminate.  S=a  ,  ...,a,  =3  with  arc  (a.,a.t1)  having 

ml  j  j+1  * 

k.  improvements  is  the  desired  path.  The  arcs  of  the  primal  upon 
J 

which  breakdowns  must  occur,  are  those  which  intersect  this  path. 
Specifically,  kj  breakdowns  must  occur  on  the  primal  arc  which 
intersects  arc  (a^.a^j)  in  this  path. 


3.  Justification  of  the  Algorithm 

The  justification  will  consist  in  showing  that  the  algorithm  is 
finite  and  that  at  termination  the  label  of  highest  rank  in  set  n  at 
3  has  as  its  first  two  components  the  mean  and  variance  of  the 
length  of  the  n-arc  path  to  the  sink  of  minimum  effective  length. 
Furthermore,  the  steps  of  the  tracing  procedure  can  be  carried  out 
and  finds  this  path. 


Lemma  1:  Let  S  =  an,  a,,  .  .  .  ,  a  with  k.  breakdowns 
-  0  1m  j 

occurring  on  arc  (a.  ,,a.)  bean  i-arc  path  to  am  such  that 

J  J 

all  of  its  arc 8  are  distinct.  Let  the  mean  and  variance  of  the  length 

of  arc  (a  .,a.)  be  u  and  'Mr?  11  respectively.  Then  the  total 
J-l  J  \J  J 


••Jb 


t 


^Note  that  p(a,  b,  0)  =  cr^{ a,  b,  0)  =  0 
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length  of  this  path  is  a  random  variable  whose  mean  is  equal  to 

m  m  2 

£  p.  and  whose  variance  is  equal  to  £  . 

j=i 1  j=i J 

Proof:  The  amounts  by  which  the  lengths  of  the  different 
arcs  canfbe  reduced  are  independent.  Hence  the  lengths  of  the 
arcs  themselves  are  independent  and  their  variances  may  be  added 
to  give  the  variance  of  the  total  length  of  the  path.  Of  course,  the 
means  may  be  added  regardless  of  whether  or  not  the  distributions 
are  independent. 


Corollary  2:  Let  S  =  a~,  .  .  .  ,  am  with  k.  breakdowns 

J 

occurring  on  arc  (a.  ,,a.)  be  a  non-cyclic  i-arc  path  to  node 

J  *  J 

a  ,  and  let  the  length  of  arc  (a.  ,,a.)  have  mean  p.  and 
m  6  j -1  y 

variance  <r?  .  Then  the  length  of  this  path  has  mean 


m 


m 


£  u.  and  variance  £  <r.  . 

M  J  j=i  3 

Proof:  Since  the  path  is  non-cyclic,  all  its  arcs  are  distinct 
and  the  corollary  follows  from  lemma  1. 

Definition:  The  pair  (p,  <r^)  is  said  to  dominate  the  pair 

—  —2 

(p,  <r  )  if  and  only  if: 


i)  M-  <  M- 

ii)  (i  +  <r  ~\Jji  -  1  <  (I  +  ^ 
with  strict  inequality  for  at  least  one  of  these.  If  path  1  is  an 


1 


i-arc  path  to  node  a  whose  length  has  mean  and  variance 

<r^  ,  and  path  2  is  an  i-arc  path  to  node  a  whose  length  has 
mean  and  variance  <r ^  then  path  1  is  said  to  dominate  path 


-  22  - 


2  if  and  only  if  the  pair  (p^.o-j)  dominates  (p^.o^)  •  An  i-arc 

path  to  node  a  that  is  not  dominated  by  any  other  i-arc  path  to 

node  a  is  said  to  be  an  undominated  i-arc  path.  Furthermore, 

a  label  (p,o-  ,  t,k)^  .  dominates  the  label  (p,  <r  ,  t,  k)^  .  if  and 

i  a,  i 

•  t"  p  *r 

only  if  the  pair  (p^  .  ,  .}  dominates  (p^  .  ,  {try  .) 

a,  i  a,  i  a,  i  a ,  i 

Note  that  the  dominance  property  is  transitive,  that  a  label 
can  be  introduced  into  a  set  only  if  no  other  label  in  that  set  domin¬ 
ates  it,  and  finally,  when  a  label  is  being  considered  for  introduction 
into  a  set,  the  labels  that  are  dropped  from  that  set  are  precisely 
those  which  are  dominated  by  the  one  being  considered. 

Lemma  3;  A  label  cannot  dominate  another  label  in  the  set. 
Proof:  Suppose  the  theorem  holds  at  one  stage  of  the  al¬ 
gorithm  and  that  node  a  and  arc  (a,b)  with  s  improvements 

■>  ; 

are  being  considered  with  respect  to  the  label  (p,  0”“,  t,  k)i  . 

0,1—8 

Let: 

^  =  Mh.  i-s  +  A*a,b)  "  H(a'b’ s) 

=  (<r2)b,i-s  +  <r2(a-b'  s> 

All  labels  among  set  i  at  node  a  whose  first  two  components 

_ 2 

are  dominated  by  (p,<r  )  are  dropped  from  this  set  and  the  con- 

dition  still  holds.  Then  the  label  (p,<r  ,t,k)^  .  =  (p,F  ,b,  s)  is 

a ,  i 

introduced  into  this  set  if  and  only  if  it  is  not  dominated  by  any  other 
label  in  this  set  and  the  property  is  still  preserved.  The  only  other 
way  a  change  in  labels  can  result  is  through  the  introduction  of 
(0,0,-,-)  or  (oe,  oo ,  - ,  -}  into  an  empty  set.  Initially,  the  con¬ 
ditions  of  the  lemma  hold  with  all  sets  empty  and  the  lemma  follows 
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from  induction. 


Lemma  4;  A  label  can  only  be  dropped  from  a  set  if  an¬ 
other  label  which  dominates  it  is  introduced  into  that  set. 

Proof:  Suppose  a  label  (p.o-  ,t,k){  is  to  be  dropped  from 

■■■  ■■  d  ,  1 

set  i  at  node  a  .  Then  there  is  a  quadruple  (p,  o-4*,  t,E)  which 

will  be  introduced  as  a  label  irto  this  set  provided  no  other  label 

would  dominate  it  and  which  has  the  property  that  it  dominates 

(|x,  <r  ,  t,k)^  Furthermore,  no  other  label  in  this  set  can  dom  in- 

a,  l 

atefitr  for  if  it  did  it  would  also  dominate  { >  t ,  .  .  Hence 

(p,  <r  ,  t,k)J  .=  (p,  a  ,7,1c)  is  introduced  into  set  i  at  node  a 
a  i  i 

proving  the  theorem. 

Corollary  5:  If  a  label  is  dropped  from  a  particular  set, 
there  will  always  be  a  label  in  that  set  which  dominates  it. 

Proof:  This  follows  immediately  from  lemma  4  and  the 
transitivity  of  the  dominance  property. 

Lemma  6:  Let  path  1  ,  S  =  aA . a  =  a  with  k.  im- 

—————  r  v  m  j 

provements  on  arc  (a._,.a  )  and  path  2  ,  S  =  aQ,  .  .  .  *am  =  a 
with  17  improvements  on  arc  (a^,a^)  by  non- cyclic  i-arc 
paths  to  node  a  Let  path  l  and  path  Z1  be  (i+r)-arc  paths  to 
node  b  formed  by  adding  arc  fa,b)  with  r  improvements  to 
path  1  and  path  2  respectively.  If  path  1  dominates  path  2  then 
path  1*  dominates  path  2' 

Proof:  Since  paths  1  and  2  are  non-cyclic,  it  follows  that 
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the  arcs  of  each  path  1,2,1',  and  2'  are  distinct.  Letting  the 
pair 8  (p(a),  <r2(a) )  .  (p(a).“ff2(a) )  ,  (p(b).  <r*fb) )  and  (p(b),  a2(b) ) 

be  the  means  and  variances  of  the  lengths  of  paths  1,  2  ,  1‘  and 
2'  ,  it  follows  from  lemma  1  that: 

p(b)  =  p(a)  +  1  {a,  b)  -  p(a,  b,  r ) 
fji(b>  -11(a)  +  l(a,bv.  -  p(a,b,rl 
<r^(b)  -  cr^(a)  +  <r^(a,b.r) 
cr^fb)  =  <r^(a)  +  cr^a,  b,  r) 

Therefore: 

p(a)  <  p>a} 

p(a)  +  i(a.b)  -  p(a.b.r)  <  p(a)  +  l(a,b)  -  p(a,b,r) 
p(b)  <  p(b) 

with  equality  only  if  \xia)  =  p(a'  and  the  first  condition  is  satisfied. 
For  the  second  condition  one  has: 

Case  1: 

p(a)  +  <r(a  )\jj}  -  1  =  p(a)  + 
p(a)  <  p(a) 
a(a)  >  o’(a) 

o2fa)  -  o^fa,  =•  fa2(a)  +  o2,a.  b,  r\]  -  [  o^fa)  +  <rZ( a,  b,  r)] 
a^(a>  -  <r^(a)  =  <r^(b)  -  <T^(b) 

<r(a)  +  "o(a)  <  ar(b)  +  a'b‘< 

<r(a)  •  a(a)  >  a,'b'  -  o\b) 
p(a)  +  [  <t( a)  -  c{a)]^Jj2  -  l  =  jl(a) 
p‘a)  +  f  <r(b)  -  afb)]^  *  l  <  p(a) 
n(b)  +  -  1  <  t(b>  +  a<b)"\|p  -  l 
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Case  2: 


p(a)  +  <rfa)  \j^  -  1  <  p(a)  +  <r( 
Define  <r  such  that 


ji(a)  +  <r  -  i  =  jJL(a)  +  <r( a;  ^  -  1 


|i(b)  + 

Also  a-(a)  <  <r 


2 - 2 — 

o-  -fa  >a 


fc'r)  \lf- 


7 

1  <  n(b)  +  a(b) 


f7' 


'hus 


<r(b) 


<  + 


ff(a,b,ri 


p(b)  +  cr(b 


■W' 


1  <  fa.(b')  -f  a(b) 


1  • 


a,  i 


2 

Lemma  7:  For  each  label  fp,<r  .t.k)1?  .  ,  the  quantities 

a ,  i 

and  (a  .  are  either  infinite  or  are  equal  to  the  mean 

fit  9  1  <k  v 


and  variance  of  an  i-arc  path  to  node  a  • 

Proof:  Suppose  at  one  stage  of  the  algorithm  that  for  each 
label  (p,  a  ,  t,  k)^  .  where  .  and  (a  ){  .  are  not  infinite, 

ct|l  at  |  1  cl  |  1 

there  is  an  i-arc  path  S  =  a^ . a^  =  a  with 


m 


c.  improve- 
J 


ments  on  arc  (a.  ,,a.)  whose  arcs  are  distinct  and  whose  length 

J  J 

has  mean  .  and  variance  (a^  and  such  that  set 
a,i  a,i 


m 


[  i  -  2  k.  ]  at  node  a  contains  a  label  whose  first  two  compon- 
j=n+l  J  n 

ents  equals  or  dominates  the  mean  and  variance  of  the  length  of  the 


path  S  =  an,  . 


,  a  with 
n 


<  breakdowns  on  arc  fa.  ,,a.) 

>  J-i  J 

Suppose  the  label  fp,a^,t,k)^  ,+f  is  introduced  through  the  ex¬ 
amination  of  arc  (a,b)  with  r  improvements  Then  there  is  a 

label  {p,a  ,t,k)^  such  that: 

a  1 1 

^b,  i+r  =  i  1  1|a  b>  *  Wa-b'r> 

■<r2)u  4i_  =  <<r2){  ;  +  <r“(a,b,r)  . 

i  i  r  d  *  i 
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Let  S  =  a q,  .  .  •  ,  =  a  with  improvements  on  arc 

(a.  |,a.)  be  a  path  which  satisfies  the  above  conditions  for 

(p,<r  ,  t,k){  .  .  Consider  the  path  S  =an . a  =  b  with  k 

a,  1  u  m+i  j 

improvements  on  arc  {a.  ,,a  )  ,  k  -  r  .  If  each  of  its  arcs 

j  m+i 

7  ~T~ 

are  distinct,  it  satisfies  the  above  conditions  for  (p,<r  ,  t.k);/  . 

r  b,  i+r 

Suppose  its  arcs  are  not  all  distinct  Then  there  is  an  n  <  m 


such  that  a  -a  and  a  .  =  a  . 

n  m  n+1  m+1 


The  mean  and  variance 


of  the  length  of  above  path  S  =■  aQ,  .  .  ,  a^+j  with  k.  improve- 

J 

m+1 

ments  on  arc  (a  ,.a  )  for  j  <  n  and  2  k.  improvements 

J  J  “  J=n+1  J 

on  arc  ^an»an+^  dominates  the  pair  tp£  i+r'^^^b  i+r^  *  A*so 
the  path  S  =  an,  .  .  .  .a  =  a  with  k.  improvements  on  arc 


n 


J 


(a^j.aj)  dominates  the  path  S  =  a^,  .  ,  am  with  k^  improve¬ 


m 


on- 


ments  on  arc  (a  ..a  )  .  Therefore,  if  2  k  =  0  there  is  a 

J  J  j=n+l  J 

^  j 

label  in  set  i  at  node  a  which  dominates  (p,  <r  ,t,k)^  .  c 

a,  i 

tradicting  lemma  3  On  the  other  hand,  suppose 

m  m 

2  k  £  0  Consider  any  label  in  set  [  i  -  2  k.  1  at  node  a 
j=n+l  J  j=n+lJ 

whose  first  two  components  dominate  the  mean  and  variance  of  the 
length  of  S  =  a^,  .  .  .  ,  a^  =  a  with  k.  improvements  on  arc 


J 


m+1 


(a.  |,a.)  This  label  and  arc  (a,  b)  with  2  k.  >  r  improve- 
J  J  j=n+l  J 

ments  will  be  examined  prior  to  the  introduction  of  (u,  <r^,  t,  k)?  . 

n  b, i+r 

as  a  label.  After  this  examination  there  must  be  a  label  among  set 

(i+r)  at  node  am+i  -  b  which  equals  or  dominates 

(p,  t,  k)£  But  this  prevents  its  introduction  as  a  label. 


LI  - 


Thus  the  arcs  of  S  =  aA, .  .  .  ,  a  . ,  =  b  arc  distinct.  The  intro- 

u  mri 

duction  of  (0,  0,  «,  -)  into  any  set  at  the  source  or  (oc,oo,  ~) 
into  any  set  preserves  the  above  properties.  Furthermore,  the 
dropping  of  a  label  preserves  these  properties  since  it  is  immedia¬ 
tely  followed  by  the  introduction  of  a  label  into  its  set  which  domin¬ 
ates  it.  Initially  these  properties  hold  with  ail  sets  of  labels  being 
empty.  The  lemma  follows  from  induction. 

Lemma  8:  The  labeling  algorithm  is  finite. 

Proof:  Since  the  first  two  components  of  a  label  are  either  in¬ 
finite  or  are  equal  to  the  mean  and  variance  of  a  path  whose  arcs 
are  distinct,  it  follows  •;ijft-this  pair  must  be  selected  from  a  finite 

set.  In  addition  the  choices  for  the  t**  and  .  are  fin  ie. 

a, i  a, i 

Thus  the  labels  themselves  are  selected  from  a  finite  set.  Further- 
more,  no  label  may  be  introduced  into  the  same  set  more  than  once 
for  if  it  is  once  dropped  there  is  another  label  in  the  set  which 
dominates  it  and  prevents  its  re-entry.  Thus  the  algorithm  is 
finite. 


Lemma  9:  At  termination  no  labels  are  infinite  at  any  node 
to  which  a  path  exists. 

Proof:  Consider  sec  i  at  node  a  .  Let  3  =  a„,  .  .  .  ,a  - 
— ~~  — ~  U  m 

a  be  any  path  to  a  It  follows  that  on  or  before  the  (m+lV^-* 
examination  of  the  nodes  /'or  set  i  in  step  5  a  finite  label  will  be 
introduced  into  the  set  i  at  node  a 
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•  2 

Corollary  10:  At  termination  .  and  (<r  .  corres- 

*  a  .  i  a  1 1 

pond  to  the  mean  and  variance  of  i>arc  paths  to  node  a  for  all  i 
and  all  j  provid  *d  a  path  to  a  exists. 

Proof:  This  follows  immediately  from  lemma  9. 


Lemma  11:  Let  S  =  aA,  .  .  .  ,  a  =  a  with  k.  improve- 
-  u  m  j 

ments  on  arc  'a  ,,a.)  be  an  undominated  i-arc  path  to  node  a 

J“ 1  J 

Then  this  path  contains  no  cycles. 

Proof:  Suppose  this  path  does  contain  a  cycle.  Let 

(a.  ^,a.)  ,  j  =  1',  .  .  .  ,m  be  the  distinct  arcs  of  this  path.  Further- 

more,  let  arc  (a  ^j.a  )  be  used  C.  times  in  this  path  and  let 

J  J  J 

its  length  have  mean  \x  and  variance  <r 7  .  Then  the  total  length 

J  J 

■  ,  m  m  2  2 

of  this  path  has  mean  Z  C.  u.  and  variance  Z  G^a  ,  C.  >  1 


J  J 


J  J  J  ~ 


all  j  .  Delete  all  cycles  from  this  r^th.  The  resulting  path  is  an 


m 


i-arc  path  to  node  a  whose  length  has  mean  Z  C  ji.  and  vari- 

j=i  j  j 


m 


r  2  2 


ance  Z  C.  <r  ,  C  <1  all  j  and  C.  =  0  for  some  j  .  Thus 
j  J  J  J  J 

the  path  S  =  aA . a  =  a  is  not  ur.dominated. 

r  0  m 


Lemma  12:  Let  S  =  aA,  .  .  .  ,  a  =  a  with  k  improve- 
-  0  m  j  r 

ments  on  arc  (a.  ,,a  )  be  an  undominated  i-ax  :  path  to  node  a 

’  J*1  J  K 


0»  •  with  k.  improvements  on  arc  (a.  j,a.) 

^  ^  J  J  J 


Then  S  =  a^,  .  .  ,  a 
is  an  undominated  (i-k  )-arc  path  to  node  a 


m  1  m-1 

Proof:  Suppose  the  above  path  is  not  undominated.  Let 

S  =  bQ,  .  '  •  =  am.^  with  It  improvements  on  arc  (b  yb  ) 

be  an  undominated  (i-k  )-arc  path  to  node  a  .  that  dominates 

m  m- 1 


•  ^9  - 


improvements  on  arc 


S  =  an . a  .  with  k 

u  m-x  j 


Then  S  =  ,  b  and  S  -  aAl...,a  , 

U  m  U  m~i 


contain  no  cycles 


and  it  follows  from  Jemma  6  that  S  s  ,  h  — ,  a  with  Ic  , 

C  m  m  ) 

improvements  on  arc  (b^  ^  b.)  and  k^  improvements  on  arc 


■  s.  on  arc  (a  T,a)  contra* 
J-l  j 


,  a  j 

is  an 

i-arc 

path 

ro  rode  a 

.  m 

m 

CL  m  •  • 

•  ^  2 

with 

k 

imo:  eve  mentis. 

0 

m 

i 

* 

dieting  the  hypothesis. 


i  2  i 

Lemma  13;  At  termination  the  set  of  pairs  «  !:  ,) 

. . . .  a,  k  a,  i 

is  the  set  of  means  and  variances  of  all  undominated  i-arc  paths 
to  node  a  provided  a  path  to  a  exists. 


Proof:  Suppose  that  after  the  algorithm  terminates  for  the. 

?  2  i 

sets  r  the  lemma  holds  for  all  se*s  of  pairs  (p/  ,  (cr  )**  )  , 

d  ;  1  &)  1 

i  <  r  .  Let  S  ~  a,  ,  .  .  .  ,  a  »  a  with  k  improvements  on  arc 


m 


(a^  ^  a^)  be  an  undomirat^d  (r+'.)-arc  path  to 


node  a  and  let 


>a  , ,a  )  be  the  last  arc  cr  this  pa*h  upon.  which  improvements 
’  p-1  p  r  r  r 

occur.  From  repeated  application  of  lemma  12  it  follows  that 

S  -  a, . a  ,  with  k  improvements  on  arc  (a  „,a  )  is  an 

G  p-1  '  y  '  '-i  v 

undominated  (r+i-k  )-arc  oath  to  nede  a  .  Let  the  mear  and 

p-  *  p -I 

2 

variance  of  the  lengths  of  these  two  paths  be  the  pairs  (u  ,  <r  } 

m  m 


and  (ut_  .J  respectively.  Then  there  exists  a  pair 

2 


“p-j/'p-i-1' 


2  * 

,  r+l-k  >1  r+L-k  }  =  Vl'V^  -  {i**"  th8reiSa 

p-I  p  p-i.  p  r  r 

label  in  set  r+l-k  at  node  a  ,  who?*  *‘is  st  two  components 

P  P~  • 

l 


are 


p-1 


And 


cr  . )  Note  that  the  path  S  --  a r  .  ...  t 

P  -  A  t 


k  improvements  on  arc.  (a  ,  a  )  le  ar  un  dominated  (r+l)-arc 

J  j-J  J 


path  to  node  ap  .  Denote  the  mean  and  variance  of  the  length  of 

this  path  by  and  <r*  There  can  be  no  pair 

j  ^  j 

,  r+l’^  ,  r+1^  that  dominates  (p  « )  since  such  a  pair 

P  p  P  P 

must  either  be  infinite  or  be  equal  to  the  mean  and  variance  of  an 


(r+l)-arc  path  to  node  .  Thus,  after  examining  set  (r+i-k  ) 
at  node  ap_j  and  arc  'ap-i»ap>  with  k^  improvements  there 
must  be  a  label  {(i>  <r2, k)^  r+.  such  that  4>r+,  =pp  and 

^a  ,  r+1  =  °p  *  Furthermore,  this  label  can  never  be  dropped. 

Repeated  application  of  this  argument  yields  the  conclusion  that 

eventually  a  label  (4,  g2 ,  t,  k)^  will  be  introduced  where 

j  2  j  2 

^a,  r+1  =  and  ^  ^a,  r+i  =  am  and  that  u  wiU  n*ver  be  dropped 

Furthermore,  any  label  in  set  r  +  i  at  node  a  whose  first  two 

components  correspond  to  the  mean  and  variance  of  a  dominated 
(r+l)-arc  path  to  node  a  will  be  dx  opped  by  the  introduction  of  a 
Uibei  whose  first  two  components  correspond  to  the  mean  and  vari¬ 


ance  of  the  length  of  a  path  that  dominates  the  above  mentioned  one 
Thus,  when  the  algorithm  terminates  for  the  sets  r+1  the  lemma 
holds  for  all  sets  of  pairs  >J  (a2)*’  )  ,  i  <  r+1  .  Let 

t  1  ex  *  1  — 

S  =  b0 . bin  =  a  be  an  undominated  O-arc  path  to  node  a  . 

An  argument  similar  to  that  above  shows  that  after,  at  most.  m+l 

checks  of  the  nodes  there  is  a  label  (4,  a2,  t,  k)J  n  such  that 
j  2  a,  P 

^a,  0  and  ^  0  arc  the  mean  and  variance  of  the  length  of 

this  path,  and  finally,  when  the  algorithm  terminates  for  the  sets 
0  the  lemma  holds  for  all  sets  of  pairs  (4^  (cr2>^  }  .  The 

”a  ,0  *  ■  a .  0 ' 
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lemma  itself  follows  from  induction. 


* 


Theorem  14:  When  the  algorithm  terminates,  the  label  of 
highest  rank  among  set  n  at  S  has  as  its  first  two  components 
the  mean  and  variance  cf  the  n-arc  path  of  minimum  effective 
length  to  15  . 

Proof:  The  n-arc  path  of  minimum  effective  length  is  undom 

inated  and  therefore  there  must  exist  a  label  (p,<r  '»t,k)g  n  at 

termination  such  that  p^j*  n  and  ^  correspond  to  the 

mean  and  variance  of  the  length  of  this  path.  Furthermore,  there 

*»  "r" 

can  bw  no  other  label  (p,  o'4-,  t.  n  of  higher  rank,  for  if  so 
one  has: 


contradicting  lemma  3. 

Of  course  the  n-arc  path  of  minimum  effective  length  to 
is  the  route  of  smallest  effective  length  possible  from  source 
to  sink  due  to  at  most  n  breakdowns.  Thus  the  algorithm  dc-es 
find  the  mean  and  variance  of  the  desired  path-  It  new  must  be 


shown  that  the  tracing  procedure  actually  finds  this  path  and  the 


arcs  upon  which  the  n  improvements  must  occur. 


Lemma  15: 

'  i.  Then  there 
such  t*  at: 


Let  (p,  o-^,  t,  k)J  .  ,  a^S  t  e  a  label  at  termina- 

is  another  label  (p,  <r  ,  t,k)£  -r-  at  termination 
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b  =  tJ 

a,  i 


i  =  i  -  kJ 


a,  1 


“i.i  *»*b.r  •  1(a’b)  ‘  Ma.b.O 


i  (ff r  +  a^(a,b,k-'  ) 

<*  •  1  D  ,  1  a  (  1 

Proof.  Consider  the  label  from  which  ( ^jl ,  cr  ,  t.k;* 

s ,  1 

was  obtained  This  label  must  satisfy  the  above  conditions,  'i.e.  , 

2m  T 

must  be  \p.<r  .t.k^y)  Furthermore,  this  label  cannot  be 

2  p 

dropped,  for  if  it  is  another  label  (p^ar  ,  t.k)^  -r-  is  introduced 

which  dominates  it.  Examining  this  label  and  arc  fa,b)  with 

k**  improvements  produces  a  label  which  dominates  and  hence 
a,  i 

drops  (p.<r^,tfk)J  . 

a ,  i 


AT 


Lemma  16:  The  tracing  procedure  is  finite. 

Proof:  The  sequence  of  labels  found  is  strictly  decreasing 
lexicographically  in  (p-1  i<  and  hence  none  can  be  repeated 

Cl;  1 

Thus  the  process  is  finite. 


Lemma  17:  Let  the  sequence  of  labels  found  by  the  tracing 
procedure  be: 


Then: 


(p,  <r2,  t,  k)*^r  .  ,  r  =  lt  .  . ,  m 

0  t  1 

r  r 


m- 1 


pi  =  p^  .  =  /  i  i  Ja  ,  a  -  p(a  .a  k  )j 
nS,n  ai»li  Lt  r  r  r  +  ‘  r 


.  2, J( 

(  /T  i  4 


r=l  , 

m-l 

A.1.  V  _2, 


’>  V.ns^  i  a  i  °  «V*r+l'kr 

r  -  ! 
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where  k  =  k  4 

r  a  ,  1 
r  r 


Proof:  It  follows  from  the  rule3  of  the  tracing  procedure 


that: 


j 

^  ^a  ,  i  "tAar+i.,i  =  1  (ar»  ar+i^  *  fA(ar  *  ar+x»  krL  r  =  I,  .  .  . ,  m-1 


r  r 


r+r  r+ i 


.  =  C 

rS, i  a  , i 
m  m  m 

■2J. 


.2, 


(ii)  (OJr  .  -  (<r2vjr+i  •*-  <r“(a  .a  ,,k  )  ,  r  =1 . m-1 

Vlr  ar+l*lr+i  r  +]  r 

(a2)Jom.  =  >Vm  •  -  0  . 

S, i  a  .  l 

m  m  m 

Summing  the  equations  (i)  give  the  first  relationship  and  summing 
the  equations  (ii)  gives  the  second. 

Theorem  18:  Let  the  sequence  cf  labels  found  by  the  tracing 

procedure  be  as  in  lemma  17.  Then  S  =  a  ,...,a,  =  5  with 

m  1 

k  =  kar  ^  improvements  on  arc  (a  ,a  .)  is  an  n-arc  path 

from  source  to  sink  whose  length  has  mean  and  variance 

b,  n 

/  2v*l 
(  “• n 

Proof:  The  quantities  \jQ  n  a-id  (o-2)^  n  equal  the  sums 
of  the  means  and  variances  respectively  of  the  arcs  in  this  path. 

All  the  arts  in  this  path  are  distinct,  for  if  not  there  would  be  a 
cycle,  and  deleting  all  cycles  would  yield  an  n-arc  path  to  15  such 
that  its  length  has  mean  strictly  less  than  jig  ^  and  variance  not 
exceeding  (<r  contradicting  theorem  14.  Thus  the  theorem 

follows  from  lemma  1. 

This  completes  the  justification  of  the  procedure  for  finding 
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the  desired  route  itself.  The  primal  arcs  upon  which  breakdowns 

must  occur  are  those  which  intersect  this  route.  In  particular  k. 

«3 

breakdowns  must  occur  upon  that  arc  of  the  primal  which  intersects 


arc  (a.,aai1)  . 

j  j+i 


4.  Normal  Case 

Suppose  it  is  known  that  the  decrease  in  the  capacities  of  the 
arcs  due  to  breakdowns  are  independently  and  normally  (or  rather 
approximately  normally)  distributed  with  known  mean  and  variance 
It  follows  that  for  any  breakdown  pattern  all  routes  through  the  dual 
are  normally  distributed. 

Let  the  length  of  a  particular  path  be  a  random  variable  L 
whose  distribution  is  normal  with  mean  p  and  variance  <r^ 

Then: 

-  £)]■  =  0 

where  $(1  -  £)  satisfies 

P  jV  <  *<l  -  0  !  j  -  1  -  0 
where  Y  is  standard  normal.  Thus: 


p|l  >  |1  +  e(i  -  0)}--  0 

Hence  p  +  <r  $(1  -  0)  plays  the  role  of  the  effective  length  of  this 
path.  Suppose  it  is  desired  to  find  the  n-arc  path  to  the  sink  of 
minimum  effective  length.  Then  the  algorithm  just  presented  in 
this  chapter  may  be  used  to  solve  this  problem  by  replacing 


*(1  -0) 
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5.  Validity  of  the  Bound 

Under  the  breakdown  pattern  found  by  the  algorithm  the 
following  relationship  exists: 

P{V[A,B]  >  F}  <  0 

where  [A,B]  is  the  cut  set  which  intersects  the  n-arc  path  of 
minimum  effective  length  through  the  dual  and  F  is  the  effective 
length  of  this  path.  However,  the  stochastic  nature  of  this  network 
indicates  that  [A,  B]  need  not  be  the  minimum  cut  with  probabil¬ 
ity  one  and  thus  one  is  led  to  believe  that  there  might  exist  an 
F'  <  F  such  that: 

P{max  flow  >  F'}  <  j3 

In  this  section,  a  stochastic  program  is  formulated  which  applies 
Tchebycheff's  extended  lemma  to  the  entire  network.  It  is  then 

shown,  surprisingly  enough,  that  such  an  F*  does  not  exist  for 
this  program. 

This  linear  programming  formulation  of  the  network  is: 

Find  x_  >  0  ,  max  F  such  that: 


X 

U. 

1  A 

b.. 

(1) 

fr* 

-5 

1  X 

.  Sj 

+  F  =  0 

X. 

j 

Yx.. 

L  Jl 

x. . 

,  ij 

=  0  .  i  i  S  ,  5 

j 

r—v 

j 

Zxjs 

-I 

xSj 

• 

II 

o 

j 

j 

The  dual  to  this  program  is: 
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Find 

(2) 


uij  ^  ®  »  vi  *  m*n  Z  such  that: 

u. .  -  v.  +  v  >0 
lJ  i  J  - 

VS  ’  v5  -  1 


=  Z 


where  b..  is  a  random  variable  Of 

"J 

theorem  of  linear  programming,  min 

values  of  the  b.  .  ^  ^  Let  u.  and 

lJ  ij 


course,  by  the  duality 

Z  -  max  F  for  any  fixed 
2. 

. .  be  the  mean  and 


variance  of  the  amount  by  which  the  capcity  of  arc  (i.j)  is  re¬ 
duced  by  the  given  breakdown  pattern  and  let  C.  be  its  original 
capacity.  Then: 


E(b. .)  =  C.  -  u.. 

VJ  iJ 

var/b. .)  =  <r^ 
iJ  lJ 

Applying  Tchebycheff* s  extended  lemma,  the  stochastic  pro 
gram  (2)  can  be  reduced  to  the  nonlinear  program  (31. 

Find  u.^  >  0,  v.,  min  Z  such  that: 


uij  +  vi  -  vj  <  0  •  (i.j)  *  <S,3> 


(3) 


-VS  +  v5  +  1  1  0 


I(C>j '  = 

’'i.j 


=  z 


I.J 


Reference  4,  pp.  128-134. 

(2) 

Reference  11,  pp.  2-6. 
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Theorem  1:  In  (3)  min  Z  <  F 


Proof:  Let 


r1 

U.  .  =  4 

lJ  lo 


1  if  i  €  A  ,  j  €  B 
otherwise 


r\  if  i  e  A 

vi  =  1 

U  if  icB 

This  satisfies  the  constraints  of  (3)  and  gives  Z  -  F  . 

It  now  remains  to  determine  whether  or  not  min  Z  =  F  .  If 
min  Z  =  F  ,  then  the  feasible  solution  to  (3)  given  in  the  proof 
of  theorem  1  is  also  an  optimal  solution.  This  happens  if  and  only 


if  there  exists  x. .  >  0  ,  F  >  0  such  that: 

ij  - 


.a) 


(4) 


G((i,0)  =  min  G(u,v)  = 


Min 
u  >  0 


z  + 


I 


(i,  j)  *  (S.'S) 


Xij(- Wvj)+  F(-vs+vsH) 


where  u  =  [  ]  ,  v  =  (v^)  and  (6,0/  is  the  .easible  solution 

to  (3)  given  in  the  proof  of  theorem  1 
This  is  satisfied  if  and  only  if: 


(5) 


1-1,  J  1  11  J 


for  all  u  >  0  at  (u,  v)  =  ((i,0)  This  happens  if  and  only  if: 


9G{u,  v) 


u. 

l 


<« 


9G(u,  v) 
dv" 


(u;  v)  =  (A*  (>) 


(u,  v)  =  (A,  0) 


(u,  v)  =  (&,  0) 


=  0  for  u  =  0 

i 


=  0  for  u.  >  0 

i 


=  0  all  i  . 


7U 


Reference  4,  pp.  471-472 
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Performing  the  differentiation  indicated  in  (6)  one  obtains: 


Ix5.  j  '  Ixj3  +  F  =  0 
j  j 


It  will  be  shown  in  the  proof  of  the  following  lemma  that  the 
relationships  (7)  can  be  satisfied  if  and  only  if  a  certain  network 
which  is  related  to  the  primal  network  has  [A,  B]  as  a  minimum 
cut.  Later  on  it  will  be  shown  that  [A.B]  is  a  minimum  cut  of 
this  network  and  therefore  (7)  can  be  satisfied  and  min  Z  :  F 

Lemma  2*  Consider  the  network  formed  by  taking  the  primal 
network  and  changing  the  capacity  of  arc  (i.j)  to 
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2 

•  (Note  that  p..  =  <r..  =  0  if  no  break- 


2 


downs  occur  on  arc  (i,  j)  ).  Then  min  Z  =  F  in  (3)  if  and  only 
if  [A,  B]  is  a  minimum  cut  of  this  network. 

Proof:  If  min  Z  =  F  in  (3)  then  the  x.^  ,  F  which  satis¬ 
fy  (7)  also  satisfy  the  constraints  of  the  network  and  produce  a 
flow  of  F  =  V[A,  B]  .  Thus  [A,B]  is  a  minimum  cut.  On  the 
other  hand,  if  [A,  B]  is  a  minimum  cut,  then,  letting  x.^  be 
the  flow  along  arc  (i,  j)  the  of  any  maximal  flow  pattern 

and  F  satisfy  (7)  . 

Corollary  3:  If  min  Z  £  F  in  (3)  then  the  network  in 
theorem  2  has  a  cut  [C,D]  such  that  V[C,D]  <  F  . 

Proof:  This  follows  immediately  from  lemma  2. 

Definition:  Let  the  effective  capacity  of  a  set  of  arcs,  A 
be  p  +  •iFi  where  p  and  cr~  are  the  iSRean  and  variance  of 
the  sum  of  the  capacities  of  the  arcs  in  A  respectively. 

Note  that  all  properties  of  effective  lengths  of  arcs  in  the 
dual  hold  for  the  effective  capacities  of  the  corresponding  arcs  in 
the  primal.  Hence  the  analogous  versions  of  the  theorems  in 
Chapter  4,  section  3,  on  effective  length  will  be  assumed  to  hold 
for  effective  capacity. 

Lemma  4:  The  number  of  breakdowns  occuring  on  cut 
[A,  B]  is  n  . 
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Proof:  Suppose  only  n  -  k  ,  k  >  0  ,  breakdowns  occur  on 


cut  [A.  B]  Increase  the  number  of  breakdowns  on  any  arc  of 
[A,  B]  by  k  .  This  results  in  a  decrease  in  the  effective  capaci¬ 
ty  of  [A,  B]  and  hence  in  the  effective  length  of  the  corresponding 
route  through  the  dual  network,  contradicting  theorems  14  and  18 
of  Chapter  4  ,  section  3 

Lemma  5:  If  min  Z  <  F  in  (3)  then  under  the  breakdown 
pattern  found  by  the  algorithm  there  is  a  cut  [C,  D]  such  that 
the  effective  capacity  of  the  arcs  of  {[C,  D]  -  [A,  B]  D  [C,  D]  }  is 
strictly  less  than  the  effective  capacity  of  the  arcs  of 
{[A,  B]  •  [A,  B)  f)  [C.D]} 

Proof:  From  corollary  3  there  exists  a  cut  [C.D]  such 
that:  , 

C1  ‘  »*1  +  <  C2  '  **2  +  '  * 

where  Cj  is  the  sum  of  the  zero-breakdown  capacities  of  the 

arcs  in  [C.Dl  ,  pj  and  are  the  mean  and  variance  of  the 

reduction  in  this  sum  due  to  the  breakdown  pattern  and  • 

1 ,  and  o’ ^  are  these  same  quantities  for  [A.B]  Let  C  , 

2 

H  ,  and  <r  be  these  quantities  for  [A,B]A[C,D]  From 
lemma  4,  p  =  p.  and  <r^  -  Thus  the  effective  capacity  of 

{[C.D]  -  [A,B]fl  [C.D]}  is  C[  -  C  and  that  of 
{[A.B]  -  [A.B]  0  [C.D]}}  is _ 

[  <c2  -  c)  -  (p2  -  nj)  ) 

and  one  has: 
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proving  the  theorem. 


Lemma  6:  The  mean  and  variance  of  the  sum  of  the  capaci¬ 
ties  of  the  arc 3  of  {jC,  D]  -  [a,  B]  [C  Dj  }  dominates  the  mean 
and  variance  of  the  sum  of  the  capacities  of  the  arcs  of 
{[A,  B]  -  [A.B]0  [C,  D]  }  . 


Proof:  Using  the  teiminology  of  theorem  5  one  has: 

-  C)  <  -  C)  -  (u?  -  (1,'j  ^  -  1 


and  since  {p2  -  -  T  <r*  -  **  l/^  -  1  >  0  it 

(C^  -  C)  <  [C^  -  C)  proving  the  theorem. 


follows  that 


Theorem  7;  In  the  convex  program  (3)  min  Z  -  T  . 

Proof:  Assume  min  Z  <  T  .  Then  there  is  a  cut  [C?  D] 
such  that  the  mean  and  variance  of  the  sum  of  the  capacities  of 
{[C,  D]  -  [A,  B]fi  [C,  D]  }  dominates  these  same  two  quantities 
for  {[A,B]  -  [A,B]f)[C,D]}  .  But  this  means  that  the  mean  and 
variance  of  V[C,  D]  dominates  that  of  V[A,  B]  ,  the  proof  of 
thie  being  identical  to  that  of  lemma  6,  chapter  4,  section  3.  This 
contradicts  the  fact  that  the  algorithm  finds  the  shortest  n-arc 
path  from  source  to  sink. 

Thus  the  solution  to  the  stochastic  program  does  not  improve 
the  bound  found  by  the  algorithm.  This  is  a  strong  argument  in 
favcr  of  the  validity  of  evaluating  a  breakdown  pattern  solely  on  the 
minimum  effective  capacity  of  ail  cut  sets. 
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6.  Practicality  of  the  Algorithm 

Since  this  algorithm  is  in  many  ways  similar  to  the  minimum 
route  algorithm^  which  is  highly  efficient,  it  would  seem  that  this 
algorithm  would  be  practical,  provided  the  number  of  labels  in  each 
set  remained  small*  While  examples  have  not  been  solved  to 
establish  whether  oi  not  the  number  of  labels  in  each  set  is  likely 
to  remain  small  there  are  indications  in  favor  of  this  happening. 

The  first  of  these  indications  is  the  assumption  that: 
p(a,  b,  i)  >  p(a,S,I) 

implies 

p(a,  b,  i)  -  <r(a,b,i)  >  p{a,T>,I)  -  (a.E.T)’^  -  1  . 

Thus,  the  difference  between  the  zero-breakdown  length  of  a  path 
and  its  effective  length  tends  to  increase  as  the  mean  decrease  in 
length  due  to  breakdowns  increases. 

The  second  of  these  indications  is  even  more  convincing.  Let 
A  and  B  be  two  paths  to  node  a  such  that  A  has  smaller 
mean  length  than  B  but  does  not  dominate  it.  Then  if  arcs  with 
breakdowns  are  added  to  the  two  paths,  the  new  path  ft  rmed  from 
A  tends  to  dominate  that  formed  from  B  This  is  summarized 
in  theorem  1. 

Theorem  1:  Consider  the  following  quantities: 

Pl  =  Pi  +  p 

?2  =  p, >  +  p 

^Reference  7,  pp. 130-134. 
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2  2  — 

Then  if  a\>a2.  ’  ^l’JjL2=^JLi‘?JL2  arK*  ^1  "  ^2  <  ^*1  "  ^ 
with  equality  if  and  only  if  or*  =  0 

Proof: 

jjl j  -  ^  -  M-£  follows  trivially 


Note  that  equality  holds  in  the  above  if  <r^  =  0  . 

The  significance  in  this  is  that  it  tends  to  keep  the  number  of 
labels  in  a  set  from  increasing  at  a  rapid  rate  as  one  gets  further 
from  the  source. 
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CHAPTER  V 


Otner  Approaches 

Two  valid  criticisms  of  the  algorithms  presented  are  that  it 
would  be  more  desirable  to  have  an  algorithm  that  worked  for  non- 
planar  networks  as  well  as  planar  networks  and  which  worked  direct¬ 
ly  with  the  primal  network.  These  problems  have  been  investigated 
and  the  purpose  of  this  section  is  to  discuss  some  of  the  difficulties 
encountered. 

A  labeling  algorithm  which  works  directly  with  the  primal 
was  developed  by  the  author  to  find  the  arc  which,  when  removed 

ill 

from  the  network,  wouid  reduce  the  maximum  flow  the  most.  'A' 
Among  other  things,  this  algorithm  required  a  list  of  all  the  arcs. 

If  this  were  to  be  extended  to  the  case  where  n  arcs  are  to  be  re¬ 
moved  (or  their  capacities  reduced),  a  list  of  all  n- tuples  of  arcs 
would  be  required.  Such  a  list  would  be  too  long  to  be  practical. 

An  undesirable  feature  of  this  problem  is  that  the  solution  for 
(n  -  1)  breakdowns  does  not  supply  useful  information  for  the  solu¬ 
tion  for  n  breakdowns  In  fact,  examples  have  been  constructed 
where  reducing  the  capacities  of  a  particular  set  of  n  arcs  reduces 
the  maximum  flow  to  zero  while  reducing  the  capacities  of  any  (n  - 1) 
arcs  in  this  set  causes  no  reduction  in  the  maximum  flow.  While 
this  produces  no  difficulties  when  working  with  routes  through  the 
dual,  it  appears  to  be  a.  insurmountable  barrier  when  working  with 
cuts  of  the  primal. 

‘^Reference  15. 
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Assuming  this  problem  cannot  be  solved  directly  using  the 
primal  network,  one  may  be  tempted  to  try  to  extend  the  notion  of 
a  dual  network  to  non-planar  networks.  A  natural  v/ay  to  do  this  is 
to  draw  the  network  on  a  sphere  and  place  an  artificial  node  wherev¬ 
er  two  arcs  that  are  net  joined  oy  a  node  intersect.  This  essentially 
replaces  each  arc  by  one  or  more  sub-arcs  in  such  a  way  as  to 
create  a  new  network  which  is  planar.  Each  sub-arc  is  assigned  a 
capacity  equal  to  that  of  the  original  arc  it  is  a  part  of.  A  cut  set 
[A,  B]  of  the  original  non-planar  network  will  be  repr  esented  in 
the  planar  network  if  and  only  if  there  is  a  cut  set  consisting  of  ex¬ 
actly  one  sub-arc  of  each  arc  in  [A,  B]  However,  some  cut 
sets  of  the  non-planar  network  may  not  be  represented  in  the  new 
planar  one.  An  example  of  this  is  the  network  of  figure  2.  The 
actual  nodes  of  the  network  are  1  ,  2  ,  3  ,  4  ,  and  5  and  the 
artificial  ones  are  a  ,  b  ,  c  .  d  ,  and  e  .  Consider  the  cut  [A,B] 
where  A  =  {1  ,  2  ,  4}  and  B  =  {3  ,  5}  Suppose  [A,B]  is 
represented  by  a  cut  [A',  B‘]  in  the  modified  network.  It  follows 
that  1  ,  2  ,  4  €  A’  and  3  ,  5  e  Bf  Since  no  sub-arc  of  (1,4) 
can  be  in  [A' ,B']  it  follows  that  c  ,  d  e  A*  and  since  no  sub¬ 
arc  of  (3,5)  is  in  [A'.B'j  it  follows  that  d  .  e  e  B‘  .  But 
this  is  impossible  since  A’/)  B'  =  0  and  therefore  [A,B]  is  not 
represented. 
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5 


FiG.  2  A  NON -PLANAR  NETWORK  WITH  CUT  SETS 
THAT  ARE  NOT  REPRESENTED  IN  THE  MODIFIED 
NETWORK.  NUMBERS  REPRESENT  THE  ACTUAL 
NODES  AND  LETTERS  THE  ARTIFICIAL  NODES. 
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Appendix  I;  Example  for  Deterministic  Case 

An  example  was  solved  using  the  deterministic  algorithm  of 
Chapter  III.  The  network  itself  is  shown  in  figure  3  and  the  re¬ 
sults  of  the  calculations  in  table  1. 
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(a,b) 

- ARC  OF  DUAL  WITH  LENGTH  a  AND 

REDUCTION  OF  LENGTH  b  IF  BREAK¬ 
DOWN  OCCURS 

-a  ARC  of  primal 

©  NODE  C  OF  DUAL 


FIG.  3  NETWORK  FOR  EXAMPLE  1 
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Table  1:  Results  of  Algorithm  for  Example  1 


I.  For  (D,t,k)a  0 


Initial  .  , 

.  y 

Iteration  1 

Iteration  2 

a 

(D.t.k) 

«*  t 

0 

(D,t,k)a>o 

a 

(D,  t,k)a,0 

S 

(0,-,-) 

s 

(0,-,-) 

S 

(0.-,-) 

1 

(oc,  -) 

1 

(l.s.0) 

1 

(1,S,0) 

2 

(oO,  “) 

2 

(l,s,0) 

2 

(1,S,0) 

3 

(oo,  -) 

3 

(2,2,0) 

3 

(2,2,0) 

4 

(oo,  -,  -) 

4 

(4,  1,0) 

4 

(4,  1,0) 

5 

(oo,  -) 

5 

(5,2,0) 

5 

(5.2,0) 

6 

(oo,  -,  -) 

6 

(4,3,0) 

6 

(4.3,0) 

S 

(oo,  -,  -) 

3 

(6,4,0) 

3 

(6,4,0) 

II.  For 

(D,  t,  k)  . 
<1  *  1 

Initial 

Iteration  1 

Iteration  2 

a 

(D.t.k) 

®  f 

1 

(D.t,k)a  i 

a 

(D,t.k)a>i 

S 

(0.-.-) 

S 

(0,-,-) 

S 

(0,-,-) 

1 

(oo.  -) 

1 

(0,s,  1) 

1 

(0,S,  1) 

2 

(oo,  -,  -) 

2 

(0,s,  1) 

2 

(0,s,  1) 

3 

(oo ,  - ,  - ) 

3 

(l.s,  1) 

3 

(l.s,  1) 

4 

(oo,  -) 

4 

(2,  1,  1) 

4 

(2,  1,  1) 

5 

(oo,-,-) 

5 

(2,2,  1) 

5 

(2,2,  1) 

6 

(oo,  -,  -) 

6 

(2,3,  1) 

6 

(2,3,  1) 

S 

(oo,  -) 

3 

(4,4,0) 

3 

(4,4,0) 

Path  is  S  ,  1,  4 , 5  with  breakdown  on  (1,4)  . 
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Appendix  II:  Example  for  Stochastic  Case 

Another  example  was  solved  using  the  stochastic  algorithm 

of  Chapter  IV-  The  network  is  shown  in  figure  4  and  the  properties 

of  the  arcs  in  table  2.  Each  arc  is  subject  to  at  most  1  breakdown. 

The  decrease  in  capacity  of  an  arc  due  to  a  breakdown  is  l{p) 

with  probability  p  and  1(1  -  p)  with  probability  (1  -  p)  • 

These  parameters  of  course  are  different  for  different  arcs.  The 

problem  was  solved  for  0  =  ^  and  the  results  summarized  in 

tables  3  and  4.  The  path  found  is  S  ,  1  ,  4 , 3  with  a  breakdown 

on  (4,3)  .  The  effective  length  of  this  path  is  -^  +  2Vre :  4,9  ■ 

—  3 

Since  a  reduction  in  the  length  of  (4,  S)  is  at  least  and  the 

u  3  .  1*1 

probability  that  this  reduction  is  net  greater  than  ^  18  3  3’ 

it  follows  that  the  true  value  of  F  is  4.5  .  Note  that  if  3  -  3 
the  effective  length  of  this  path  and  the  true  value  of  F  for  this 
path  (but  not  necessarily  for  the  entire  network)  are  both  4.5 
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Table  2:  Data  for  Network  of  Figure  4 


Arc 

£ 

l(l-p) 

l(a,b) 

H  (a,  b,  1) 

ff2(a.b,l) 

(S,  1) 

1/2 

1 

3/4 

1 

7/8 

1/64 

(S,2) 

1/4 

5/8 

3/4 

1 

23/32 

13/4096 

(S,  3) 

1/3 

4/3 

11/6 

3 

5/3 

1/18 

(1.2) 

1/4 

5 

4 

6 

17/4 

3/16 

(1,4) 

1/5 

1 

7/4 

3 

8/5 

9/100 

(2,3) 

1/2 

1/2 

3/4 

1 

5/8 

1/64 

(2,5) 

1/5 

2 

3 

4 

14/5 

4/25 

(3.6) 

1/6 

3/2 

1 

2 

13/12 

5/144 

14,5) 

1/3 

1 

3/2 

2 

4/3 

1/18 

(4.5) 

1/3 

3/2 

2 

2 

11/6 

1/18 

(5,  6) 

1/4 

4 

3 

5 

13/4 

3/16 

(6.?) 

1/2 

7/4 

9/4 

3 

2 

1/16 

P  *■=  1/5 
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Table  3:  Results  of  Algorithm  for  Example  2  -  i  =  0 


Initial  Iteration  1 


a 

<r2.t.k)^0 

a 

(p,  <r2.  t,k) 

S 

(0.0,-,-) 

s 

(0,0,-,-) 

1 

(eo,  00,  .,  .) 

1 

(1,0,S,0) 

2 

(oc,  00,  -,  -) 

2 

(1,0,S,0) 

3 

(oo,  00,  -,  -) 

3 

(2, 0,2,0) 

4 

(oo,  00,  -  ,  .) 

4 

(4,0,  1,0) 

5 

(oo,  00,  .,  -) 

5 

(5, 0,2,0) 

6 

(oo ,00,  -,  -) 

6 

(4,  0,3,0) 

3 

(oo,  00,  -,  -) 

3 

(6, 0,4,0) 

Iteration  2 
a  (p,  o-  »t,k)^  0 
S  (0,0,-,-) 

1  (1.0,  S.O) 

2  (1,0,  S,0) 

3  (2, 0,2,0) 

4  (4,0,  1,0) 

5  (5, 0f  2,  0) 

6  (4, 0,3,0) 

3  (6, 0,4,0) 
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Table  4:  Results  of  Algorithm  for  Example  2  -  i  =  1 


a 

(p,<r2,t,k)^  x 

a  (|x,  <r  ,  t.k)^  j 

S 

(0,0,-,-) 

S  (0,0,-,-) 

1 

(oo,  oc,  -) 

1  <1/8,  1/64. S,  1) 

2 

(oo,  00,  -  ,  -) 

2  (9/32,  13/4096,  S,  1) 

3 

(oo,  00,  -) 

3  (41/32,13/4096,2,0) 

4 

(oo,  00,  -,  -) 

4  (12/5,9/100,1.1) 

5 

(oo,  00,  -,  -) 

5(  (11/5,4/25,2,1) 

6 

(oo,  00,  -) 

6  (35/12,5/144,3,1) 

5 

(oo,  00,  -,  -) 

S  (25/6,1/18,4,1) 

a 

•y 

(p,<r  ,  t.k)^  j 

s 

(0,  0,  - r t) 

1 

(l/8,cl/64,,S,  1) 

2 

(9/32, 13/4096, S,  1) 

3 

(41/32,  13/4096, S,  1) 

4 

(12/5,9/100,1,  1) 

5 

(11/5,4/25,2,  1) 

6 

(35/12,5/144,3,  1) 

3 

(25/6,  1/18,4,1) 

Path  is  S 

,  1,4,3 

with  breakdown  on  (4, *5)  . 

Effective  length  of  path  is  approximately  4.9  . 
True  value  is  4.  5  . 
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LIST  OF  SYMBOLS 


(a,b) 

l(a.b) 

d{a,b) 


Axe  joining  a  and  b 

Length  of  arc  joining  a  and  b 

Deterministic  decrease  in  length  of  (a,b)  resulting 
from  a  breakdown 


(D,  t,  k)  .A  label  for  the  deterministic  case 
a,  1 


D  . 
a,  1 


k 

a,  1 
L 

a,  i 


First  component  of  (D,  t,  k)  . 

1 

Second  component  of  (D,t,k)  . 

3Ly  1 

Third  component  of  (D,  t,  k)  . 

Length  of  the  shortest  i-arc  path  from  the  source  to 
node  a 


p(a,b,  i;  Mean  of  the  decrease  in  capacity  of  (a,  b)  due  to  i 

breakdowns 

cr^(afb,ij  Variance  of  the  decrease  in  capacity  of  (a,b)  due  to 

breakdowns 


(p.  o-  ,t,k^  * 

-4.1  V 

s< 

4.  t  T 

.  F 

a,  i 

[A,  B]  C 

V[A,  B]  V 


A  labei  for  the  stochastic  case 

2  i 

First  component  of  {p,<r  ,  t,k)J  . 

a,  x 

2  * 

Second  component  of  {p,  <r  ft,k}{  . 

a,  i 

Third  component  of  (p,  <r^,t,  k)^  . 

a,  i 

Fourth  component  of  (p,  c  ,  t,k){  . 
Cut  set  [A,  B] 

Value  of  cut  set  [A,  B] 


REFERENCES 


1.  Boldyreff ,  A-  W.  ,  Determination  of  the  Maximal  Steady- 
State  Flow  of  Traffic  Through  a  Railroad  Network,  'The  RAND 
Corporation,  ft.M-1532  (ASTlA  No.  88712),  August  5,  1955. 

2.  Brossman,  M.  W.  ,  R.  G.  Busacker,  et  al,  Computer  Assist* 

ed  Strategic  Logistic  Planning  Transportation  Phase,  Opera- 
n"6nS“ITe  search  Office",  Jul”176i: - 

3.  Busacker,  R.  G*  ,  and  P.  J.  Gowen,  A  Procedure  for  Deter¬ 
mining  a  Family  of  Minimum  Cost  Network  Flow  Patterns, 
Operations  Research  Office  0RO-TP-l5,  November,  I960. 

4.  Dantzig,  G.  B.  ,  Linear  Programming  and  Extensions, 
Princeton  University  Press,  Princeton,  New  Jersey,  1963. 

5.  Dantzig,  G.  B.  ,  and  D.  R.  Fulkerson,  Notes  on  Linear 

Prog  ramming --Part  XXI:  On  the  Max-Flow  Min-Cut  Theorem 
of  Networks,  The  RaND  Corporation,  RM-i4l8-l  (AsTIA  No- 
86705),  January  1,  1955. 

6.  Dantzig,  G.  B.  ,  and  D.  L.  Johnson,  Maximum  Payloads  per 
Unit  Time  Delivered  Through  an  Air  Network,  Boeing  Scien¬ 
tific  Research  Laboratories,  Df- 82-02 6o,  June,  1963. 

7.  Ford,  L.  R.  ,  Jr.  ,  and  D.  R.  Fulkerson,  Flows  in  Networks, 
The  RAND  Corporation,  R-375,  December  20,  I960. 

8.  ."Maximal  Flow  Through  a  Network,"  Can.  Jour. 
Rath.  .  TJ  (1956),  399-404. 


9-  Frechet,  M. ,  Recherches  Theoriques  Modernes  sur  la 
Theorie  des  Probabilities,  Paris,  19^7. 

10.  Fulkerson,  D.  R.  ,  An  Out- of- Kilter  Method  for  Minimal  Cost 

Flow  Problems,  The  RArdD  Corporation,  P-1825,  April  5, 
TTOT - 

11.  Sinha,  S.  M- ,  Stochastic  Programming,  ORC  63-22  (RR), 
Operations  Research  Center,  University  of  California, 
Berkeley,  August  19,  1963. 

12.  Whitney,  H.  ,"Non-Separable  and  Planar  Graphs, "Trans. 
Amer.  Math.  Soc.  ,  34  (1932),  339*362. 

13.  L.,  '  ,  "Planar  Graphs,"  Fundamenta  Mathematicae 

(1933),  73-84.  - 


-  57  - 


14.  Wollmer,  R.  D- ,  Removing  Arcs  from  a  Network, 

ORC  64-5  (RR),  Operations  Research  Center,  University 
of  California,  Berkeley,  March  13,  1964. 

15.  ,  Some  Methods  for  Determining  the  Most  Vital 
Link  in  a  Railway  Network,  'fhe  RaND  Corporation, 
&hl-33Zl-£SA,  April,  19b3. 


-  58  - 


% 


